Unsupervised Segmentation of 5D Hyperpolarized Carbon-13 MRI Data Using a Fuzzy Markov Random Field Model

Hyperpolarized MRI with 13C-labelled compounds is an emerging clinical technique allowing in vivo metabolic processes to be characterized non-invasively. Accurate quantification of 13C data, both for clinical and research purposes, typically relies on the use of region-of-interest analysis to detect and compare regions of altered metabolism. However, it is not clear how this should be determined from the five-dimensional data produced and most standard methodologies are unable to exploit the multidimensional nature of the data. Here we propose a solution to the novel problem of 13C image segmentation using a hybrid Markov random field model with continuous fuzzy logic. The algorithm fully utilizes the multi-dimensional data format in order to classify each voxel into one of six distinct classes based on its metabolic characteristics. Bayesian priors fully incorporate spatial, temporal and ratiometric contextual information whilst image contrast from multiple spectral dimensions are considered concurrently by using an analogy from color image segmentation. Performance of the algorithm is demonstrated on in silico data where the superiority of the approach over a reference thresholding method is consistently observed. Application to in vivo animal data from a pre-clinical subcutaneous tumor model illustrates the ability of the MRF algorithm to successfully detect tumor location whilst avoiding image artefacts. This work has the potential to assist the analysis of human hyperpolarized 13C data in the future.CJD is jointly funded by the National Institute for Health Research (NIHR), Cambridge Biomedical Research Centre and GlaxoSmithKline (GSK). The authors acknowledge further research support from Cancer Research UK (C19212/A911376, C19212/A16628), the Cancer Research UK/Engineering and Physical Sciences Research Council Imaging Centre in Cambridge and Manchester, the CRUK Cambridge Centre and Cambridge Experimental Cancer Medicine Centre

Clinical significance of prostate (18)F-labelled fluorodeoxyglucose uptake on positron emission tomography/computed tomography: A five-year review

AIM: To determine the significance and need for investigation of incidental prostatic uptake in men undergoing (18)F-labelled fluorodeoxyglucose ((18)F-FDG) positron emission tomography/computed tomography (PET/CT) for other indications. METHODS: Hospital databases were searched over a 5-year period for patients undergoing both PET/CT and prostate magnetic resonance imaging (MRI). For the initial analysis, the prostate was divided into six sectors and suspicious or malignant sectors were identified using MRI and histopathology reports respectively. Maximum and mean (18)F-FDG standardised uptake values were measured in each sector by an investigator blinded to the MRI and histopathology findings. Two age-matched controls were selected per case. Results were analysed using a paired t-test and one-way ANOVA. For the second analysis, PET/CT reports were searched for prostatic uptake reported incidentally and these patients were followed up. RESULTS: Over a 5-year period, 15 patients underwent both PET/CT and MRI and had biopsy-proven prostate cancer. Malignant prostatic sectors had a trend to higher (18)F-FDG uptake than benign sectors, however this was neither clinically nor statistically significant (3.13 ± 0.58 vs 2.86 ± 0.68, P > 0.05). (18)F-FDG uptake showed no correlation with the presence or histopathological grade of tumour. (18)F-FDG uptake in cases with prostate cancer was comparable to that from age-matched controls. Forty-six (1.6%) of 2846 PET/CTs over a 5-year period reported incidental prostatic uptake. Of these, 18 (0.6%) were investigated by PSA, 9 (0.3%) were referred to urology, with 3 (0.1%) undergoing MRI and/or biopsy. No cases of prostate cancer were diagnosed in patients with incidental (18)F-FDG uptake in our institute over a 5-year period. CONCLUSION: (18)F-FDG uptake overlaps significantly between malignant and benign prostatic conditions. Subsequent patient management was not affected by the reporting of incidental focal prostatic uptake in this cohort

The Prostate Health Index adds predictive value to multi-parametric MRI in detecting significant prostate cancers in a repeat biopsy population

Both multi-parametric MRI (mpMRI) and the Prostate Health Index (PHI) have shown promise in predicting a positive biopsy in men with suspected prostate cancer. Here we investigated the value of combining both tests in men requiring a repeat biopsy. PHI scores were measured in men undergoing re-biopsy with an mpMRI image-guided transperineal approach (n = 279, 94 with negative mpMRIs). The PHI was assessed for ability to add value to mpMRI in predicting all or only significant cancers (Gleason ≥7). In this study adding PHI to mpMRI improved overall and significant cancer prediction (AUC 0.71 and 0.75) compared to mpMRI + PSA alone (AUC 0.64 and 0.69 respectively). At a threshold of ≥35, PHI + mpMRI demonstrated a NPV of 0.97 for excluding significant tumours. In mpMRI negative men, the PHI again improved prediction of significant cancers; AUC 0.76 vs 0.63 (mpMRI + PSA). Using a PHI≥35, only 1/21 significant cancers was missed and 31/73 (42%) men potentially spared a re-biopsy (NPV of 0.97, sensitivity 0.95). Decision curve analysis demonstrated clinically relevant utility of the PHI across threshold probabilities of 5-30%. In summary, the PHI adds predictive performance to image-guided detection of clinically significant cancers and has particular value in determining re-biopsy need in men with a negative mpMRI

The longitudinal effect of ejaculation on seminal vesicle fluid volume and whole-prostate ADC as measured on prostate MRI.

OBJECTIVE: To prospectively investigate the longitudinal effect of ejaculatory abstinence on MRI-measured seminal vesicle (SV) volume and whole-prostate ADC over consecutive days. METHODS: 15 healthy male volunteers (mean 35.9 years, range 27-53) underwent 3-T MRI at baseline and 1, 2 and 3 days post-ejaculation. Prostate and SV volumes were derived by volume segmentation and whole-gland apparent diffusion coefficient (ADC) values calculated. A mixed-effects linear regression compared ADC values and prostate/seminal vesicle volumes in each volunteer between studies in a pairwise manner. RESULTS: All subjects completed the four MRIs. Mean prostate volume was 22.45 cm³ (range 13.04-31.21 cm³), with no change between the four studies (p = 0.89-0.99). 13/15 subjects showed SV volume reduction from baseline to day 1, with group-mean decreasing from 6.45 to 4.80 cm³ (-25.6%, p < 0.001), and a significant reduction from baseline to day 2 (-18.1%, p = 0.002). There was a significant volume increase from both day 1 (+21.3%, p = 0.006) and day 2 (+10.2%, p = 0.022) to day 3 post-ejaculation. There was a significant reduction in ADC from 1.105 at baseline to 1.056 × 10¯³ mm²/s at day 1 (mean -4.3%, p = 0.009). CONCLUSION: The longitudinal effect of ejaculation on SV volume was demonstrated. Significant reductions in SV volume and whole-gland ADC were observed post-ejaculation, supporting a 3-day period of abstinence before prostate MRI. KEY POINTS: • Seminal vesicle volume significantly reduced 24 h post-ejaculation remaining reduced at day 2 • Seminal vesicle fluid volume significantly increased from day 1 to day 3 post-ejaculation • There was a significant reduction in whole-gland prostate ADC values day 1 post-ejaculation • 3-day abstinence from ejaculation is required to ensure maximal seminal vesicle distension.The authors acknowledge grant support from the Royal College of Radiologists UK and research support from Cancer Research UK, National Institute of Health Research Cambridge Biomedical Research Centre, Cancer Research UK and the Engineering and Physical Sciences Research Council Imaging Centre in Cambridge and Manchester and the Cambridge Experimental Cancer Medicine Centre

Cardiovascular Applications of Hyperpolarized MRI

Many applications of MRI are limited by an inherently low sensitivity. Previous attempts to overcome this insensitivity have focused on the use of MRI systems with stronger magnetic fields. However, the gains that can be achieved in this way are relatively small and increasing the magnetic field invariably leads to greater technical challenges. More recently, the development of a range of techniques, which can be gathered under the umbrella term of “hyperpolarization,” has offered potential solutions to the low sensitivity. Hyperpolarization techniques have been demonstrated to temporarily increase the signal available in an MRI experiment by as much as 100,000-fold. This article outlines the main hyperpolarization techniques that have been proposed and explains how they can increase MRI signals. With particular emphasis on the emerging technique of dynamic nuclear polarization, the existing preclinical cardiovascular applications are reviewed and the potential for clinical translation is discussed

Removing rician bias in diffusional kurtosis of the prostate using real-data reconstruction.

Purpose To compare prostate diffusional kurtosis imaging (DKI) metrics generated using phase-corrected real data with those generated using magnitude data with and without noise compensation (NC).Methods Diffusion-weighted images were acquired at 3T in 16 prostate cancer patients, measuring 6 b-values (0-1500 s/mm2 ), each acquired with 6 signal averages along 3 diffusion directions, with noise-only images acquired to allow NC. In addition to conventional magnitude averaging, phase-corrected real data were averaged in an attempt to reduce rician noise-bias, with a range of phase-correction low-pass filter (LPF) sizes (8-128 pixels) tested. Each method was also tested using simulations. Pixelwise maps of apparent diffusion (D) and apparent kurtosis (K) were calculated for magnitude data with and without NC and phase-corrected real data. Average values were compared in tumor, normal transition zone (NTZ), and normal peripheral zone (NPZ).Results Simulations indicated LPF size can strongly affect K metrics, where 64-pixel LPFs produced accurate metrics. Relative to metrics estimated from magnitude data without NC, median NC K were lower (P < 0.0001) by 6/11/8% in tumor/NPZ/NTZ, 64-LPF real-data K were lower (P < 0.0001) by 4/10/7%, respectively.Conclusion Compared with magnitude data with NC, phase-corrected real data can produce similar K, although the choice of phase-correction LPF should be chosen carefully

Prospective study evaluating the relative sensitivity of 18F-NaF PET/CT for detecting skeletal metastases from renal cell carcinoma in comparison to multidetector CT and 99mTc-MDP bone scintigraphy, using an adaptive trial design.

BACKGROUND: The detection of occult bone metastases is a key factor in determining the management of patients with renal cell carcinoma (RCC), especially when curative surgery is considered. This prospective study assessed the sensitivity of (18)F-labelled sodium fluoride in conjunction with positron emission tomography/computed tomography ((18)F-NaF PET/CT) for detecting RCC bone metastases, compared with conventional imaging by bone scintigraphy or CT. PATIENTS AND METHODS: An adaptive two-stage trial design was utilized, which was stopped after the first stage due to statistical efficacy. Ten patients with stage IV RCC and bone metastases were imaged with (18)F-NaF PET/CT and (99m)Tc-labelled methylene diphosphonate ((99m)Tc-MDP) bone scintigraphy including pelvic single photon emission computed tomography (SPECT). Images were reported independently by experienced radiologists and nuclear medicine physicians using a 5-point scoring system. RESULTS: Seventy-seven lesions were diagnosed as malignant: 100% were identified by (18)F-NaF PET/CT, 46% by CT and 29% by bone scintigraphy/SPECT. Standard-of-care imaging with CT and bone scintigraphy identified 65% of the metastases reported by (18)F-NaF PET/CT. On an individual patient basis, (18)F-NaF PET/CT detected more RCC metastases than (99m)Tc-MDP bone scintigraphy/SPECT or CT alone (P = 0.007). The metabolic volumes, mean and maximum standardized uptake values (SUV mean and SUV max) of the malignant lesions were significantly greater than those of the benign lesions (P < 0.001). CONCLUSIONS: (18)F-NaF PET/CT is significantly more sensitive at detecting RCC skeletal metastases than conventional bone scintigraphy or CT. The detection of occult bone metastases could greatly alter patient management, particularly in the context when standard-of-care imaging is negative for skeletal metastases.This work was supported by Cancer Research UK [grant number C19212/A16628]. The authors also received research support from the National Institute of Health Research Cambridge Biomedical Research Centre, Engineering and Physical Sciences Research Council Imaging Centre in Cambridge and Manchester, and the Cambridge Experimental Cancer Medicine Centre. The research has also been partly funded by a generous donation from the family and friends of a patient.This is the final version of the article. It first appeared from Oxford University Press via http://dx.doi.org/10.1093/annonc/mdv28

Feasibility of metabolic imaging of hyperpolarized 13C-pyruvate in human breast cancer

Introduction Imaging of the breast with hyperpolarized 13C yields new challenges compared to imaging the prostate [1]. E.g. large anteroposterior B0 gradients [2] require correction and the anatomy and patient positioning need a new, highly optimized RF coil array for achieving sufficient SNR/spatial resolution. As a first step, we have investigated single-breast imaging in the coronal plane. Methods A BRCA gene carrier with a 38-mm diameter grade 3 triple-negative invasive ductal carcinoma was studied on a 3T MRI (GE Healthcare) using a prototype 8-channel 13C breast coil (Rapid Biomedical), containing 2 transmit/receive coils and 6 receive-only covering both breasts in a prone position. 1H imaging was performed with the body coil. Following injection of 40ml of 250mM 13C-pyruvate, polarized to c. 25%, a 1-minute time series of spirals with IDEAL encoding (3) was collected (flip angle 10°, TR=260ms, 8-step cycle, time resolution 2.08s, 3 x 3-cm thick slices, 3mm gap, 40-pt spiral, 24cm coronal FOV, real pixel size 12 x 12 x 30mm). IDEAL reconstruction of images was optimized separately for each slice to enable independent frequency offsets to be applied. Kinetic modelling was performed in MATLAB, with automated tumour segmentation. Results Tumour pixels were identified by the segmentation algorithm only in the tumour-containing slice 2, and the average estimated flux from pyruvate to lactate kPL within this ROI was 0.022 s-1 (Fig. 1). The frequency shift of pyruvate relative to slice 2 was +6 Hz in slice 3 and -34 Hz in slice 1, confirming a sharp gradient in B0 approaching the nipple, which was corrected by optimizing slices separately (Fig 2). Images of lactate and pyruvate summed over the time course (Fig 3) showed strong signal of both metabolites over the tumour in slice 2, lower pyruvate in the slice toward the chest wall, and no consistent signal in slice 1. Conclusion This first-in-Europe study in breast cancer established the feasibility of obtaining metabolite images with high temporal and moderate spatial resolution in humans in vivo following administration of hyperpolarized 13C-pyruvate. Coronal image orientation allowed application of significant corrections for a known limitation, the anteroposterior B0 gradient, as well as a small FOV to improve spatial resolution. Kinetic rate constants within the tumour were found to be consistent with previous reports in human prostate cancer (1). References 1) Nelson SJ et al. Sci Transl Med 5, 198ra108 (2013). 2) Maril N, et al. Magn. Reson. Med. 2005; 54:1139-1145. 3) Wiesinger F, et al. Magn Reson Med 2012; 68:8-16

A computational model of open-irrigated radiofrequency catheter ablation accounting for mechanical properties of the cardiac tissue

Radiofrequency catheter ablation (RFCA) is an effective treatment for cardiac arrhythmias. Although generally safe, it is not completely exempt from the risk of complications. The great flexibility of computational models can be a major asset in optimizing interventional strategies, if they can produce sufficiently precise estimations of the generated lesion for a given ablation protocol. This requires an accurate description of the catheter tip and the cardiac tissue. In particular, the deformation of the tissue under the catheter pressure during the ablation is an important aspect that is overlooked in the existing literature, that resorts to a sharp insertion of the catheter into an undeformed geometry. As the lesion size depends on the power dissipated in the tissue, and the latter depends on the percentage of the electrode surface in contact with the tissue itself, the sharp insertion geometry has the tendency to overestimate the lesion obtained, especially when a larger force is applied to the catheter. In this paper we introduce a full 3D computational model that takes into account the tissue elasticity, and is able to capture the tissue deformation and realistic power dissipation in the tissue. Numerical results in FEniCS-HPC are provided to validate the model against experimental data, and to compare the lesions obtained with the new model and with the classical ones featuring a sharp electrode insertion in the tissue.La Caixa 2016 PhD grant to M. Leoni, and Abbott non-conditional grant to J.M. Guerra Ramo

Reproducibility of acute pulmonary vein isolation guided by the ablation index.

BACKGROUND: Atrial fibrillation (AF) ablation outcome is still operator dependent. Ablation Index (AI) is a new lesion quality marker that has been demonstrated to allow acute durable pulmonary vein (PV) isolation followed by a high single-procedure arrhythmia-free survival. This prospective, multicenter study was designed to evaluate the reproducibility of acute PV isolation guided by the AI. METHODS: A total of 490 consecutive patients with paroxysmal (80.4%) and persistent AF underwent first time PV encircling and were divided in four study groups according to operator preference in choosing the ablation catheter (a contact force [ST] or contact force surround flow [STSF] catheter) and the AI setting (330 at posterior and 450 at anterior wall or 380 at posterior and 500 at anterior wall). Radiofrequency was delivered targeting interlesion distance ≤6 mm. RESULTS: The rate of first-pass PV isolation (ST330 90 ± 16%, ST380 87 ± 19%, STSF330 90 ± 17%, STSF380 91 ± 15%, P = .585) was similar among the four study groups, whereas procedure (ST330 129 ± 44 minutes, ST380 144 ± 44 minutes, STSF330 120 ± 72 minutes, STSF380 125 ± 73 minutes, P < .001) and fluoroscopy time (ST330 542 ± 285 seconds, ST380 540 ± 416 seconds, STSF330 257 ± 356 seconds, STSF380 379 ± 454 seconds, P < 0.001) significantly differed. The difference in the rate of first-pass isolation was not statistical different (P = .06) among the 12 operators that performed at least 15 procedures. CONCLUSIONS: An ablation protocol respecting strict criteria for contiguity and quality lesion results in high and comparable rate of acute PV isolation among operator performing ablation with different catheters, AI settings, procedure, and fluoroscopy times